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Abstract 

This paper surveys recent work on Lie algebras of differential operators and their 
application to the construction of quasi-exactly solvable Schrodinger operators. 

1 Introduction 

Lie-algebraic and Lie group theoretic methods have played a significant role in the de- 
velopment of quantum mechanics since its inception. In the classical applications, the 
Lie group appears as a symmetry group of the Hamiltonian operator, and the associated 
representation theory provides an algebraic means for computing the spectrum. Of par- 
ticular importance are the exactly solvable problems, such as the harmonic oscillator or 
the hydrogen atom, whose point spectrum can be completely determined using purely 
algebraic methods. The fundamental concept of a "spectrum generating algebra" was 
introduced by Arima and Iachello, 0, j7|, to study nuclear physics problems, and sub- 
sequently, by Iachello, Alhassid, Giirsey, Levine, Wu and their collaborators, was also 
successfully applied to molecular dynamics and spectroscopy, |23], Ppl, and scattering 



theory, ||, [|J, ||. The Schrodinger operators amenable to the algebraic approach as- 
sume a Lie- algebraic form, meaning that they belong to the universal enveloping algebra 
of the spectrum generating algebra. Lie-algebraic operators reappeared in the discovery 
of Turbiner, Shifman, Ushveridze, and their collaborators, fl3l| , p3| , H34| , Hi) °f a new 
class of physically significant spectral problems, which they named quasi-exactly solvable, 
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having the property that a (finite) subset of the point spectrum can be determined using 
purely algebraic methods. This is an immediate consequence of the additional require- 
ment that the hidden symmetry algebra preserve a finite-dimensional representation space 
consisting of smooth wave functions. In this case, the Hamiltonian restricts to a linear 
transformation on the representation space, and hence the associated eigenvalues can be 
computed by purely algebraic methods. Connections with conformal field theory, |21 



|2T| , [R2L pi] , and the theory of orthogonal polynomials, |J7||, lend additional 

impetus for the study of such problems. 

In higher dimensions, much less is known than in the one- dimensional case; in fact, 
only a few special examples of quasi-exactly solvable problems in two dimensions have 
appeared in the literature to date, f33fl , fl8f . Complete lists of finite-dimensional Lie alge- 
bras of differential operators are known in two (complex) dimensions; there are essentially 
24 different classes, some depending on parameters. The quasi-exactly solvable condition 
imposes a remarkable quantization constraint on the cohomology parameters classifying 
these Lie algebras. This phenomenon of the "quantization of cohomology" has recently 
been given an algebro-geometric interpretation, (JT2J. Any of the resulting quasi-exactly 
solvable Lie algebras of differential operators can be used to construct new examples 
of two-dimensional quasi-exactly solvable spectral problems. An additional complica- 
tion is that, in higher dimensions, not every elliptic second-order differential operator is 
equivalent to a Schrodinger operator (i.e., minus Laplacian plus potential), so not every 
Lie-algebraic operator can be assigned an immediate physical meaning. The resulting 
"closure conditions" are quite complicated to solve, and so the problem of completely 
classifying quasi-exactly solvable Schrodinger operators in two dimensions appears to be 



too difficult to solve in full generality. A variety of interesting examples are given in [JT8 
and we present a few particular cases of interest here. 

The above ideas, originally introduced for scalar Hamiltonians describing spinless par- 
ticles, can be generalized to include spin. The first step in this direction was taken by 



Shifman and Turbiner, | 33| , using the fact that a Hamiltonian for a spin 1/2 particle in 
d spatial dimensions can be constructed from a Lie superalgebra of first order differen- 
tial operators in d ordinary (commuting) variables and one Grassmann (anticommuting) 
variable. Alternatively, 0,2x2 matrices (or N x N matrices for particles of arbitrary 
spin, 0) can be used to represent the Grassmann variable. In contrast with the scalar 
case, very few examples of matrix quasi-exactly solvable Schrodinger operators have been 



found thus far, P5| |. There are important conceptual reasons for this fact. First of all, 
in the matrix case Lie superalgebras of matrix differential operators come naturally into 
play, while in the scalar case only Lie algebras need be considered. Secondly, as we shall 
see in Section |5], every scalar second order differential operator in one dimension can be 
transformed into a Schrodinger operator of the form —d^ + V(x) by a suitable change 
of the independent variable x and a local rescaling of the wave function. For matrix 
differential operators, on the other hand, the analogue of this result — V(x) being now 
a Hermitian matrix of smooth functions — is no longer true unless the operator satisfies 
quite stringent conditions. 

In the last section of this paper, we will discuss the characterization of the class of 
quasi-exactly solvable matrix differential operators preserving a finite-dimensional space 
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of wave functions with polynomial components (cf. [|TJ, 0). For the important 

particular case of spin 1/2 particles, we will give necessary and sufficient conditions for 
a quasi-exactly solvable operator to be equivalent to a non-trivial Schrodinger operator. 
A suitable simplification of these conditions will then be used to construct new examples 
of multi-parameter quasi-exactly solvable spin 1/2 Hamiltonians in one dimension. 

2 Quasi-Exactly Solvable Schrodinger Operators 

Let M denote an open subset of Euclidean space M. d with coordinates x = (x 1 , . . . ,x d ). 
The time-independent Schrodinger equation for a differential operator 7i is the eigenvalue 
problem 

n-i) = \i). (i) 

In the quantum mechanical interpretation, 7i is a (self-adjoint) second-order differen- 
tial operator, which plays the role of the quantum "Hamiltonian" of the system. A 
nonzero wave function %p(x) is called em normalizable if it is square integrable, i.e., lies 
in the Hilbert space L 2 (M d ), and so represents a physical bound state of the quantum 
mechanical system, the corresponding eigenvalue determining the associated energy level. 
Complete explicit lists of eigenvalues and eigenfunctions are known for only a handful of 
classical "exactly solvable" operators, such as the harmonic oscillator. For the vast ma- 
jority of quantum mechanical problems, the spectrum can, at best, only be approximated 
by numerical computation. The quasi-exactly solvable systems occupy an important in- 
termediate station, in that a finite part of the spectrum can be computed by purely 
algebraic means. 

To describe the general form of a quasi-exactly solvable operator, we begin with a 
finite-dimensional Lie algebra g spanned by r linearly independent first-order differential 
operators 

r = Y^C(x)^- + V a (x), a = l,...,r, (2) 

whose coefficients ^ ai ^r] a are smooth functions of x. Note that each differential operator 
is a sum, J a = v a + rj a , of a vector field v a = ^^ ai d / dx l (which may be zero) and a 
multiplication operator 7] a . 

A differential operator is said to be Lie-algebraic if it lies in the universal enveloping 
algebra U(q) of the Lie algebra q, meaning that it can be expressed as a polynomial in 
the operators J a . In particular, a second-order differential operator T is Lie-algebraic if 
it can be written as a quadratic combination 

- T = CabJ a J b + E C * ja + C 0, (3) 

a,b a, 

for certain constants c a b, c a , Cq. (The minus sign in front of the Hamiltonian is taken 
for later convenience.) Note that the commutator [J C ,T] of the Hamiltonian with any 
generator of g, while still of the same Lie-algebraic form, is not required to vanish. There- 
fore, the hidden symmetry algebra q is not a symmetry algebra in the traditional sense. 



LIE ALGEBRAS AND PARTIAL INTEGRABILITY 



4 



Lie-algebraic operators appeared in the early work of Iachello, Levine, Alhassid, Giirsey 
and collaborators in the algebraic approach to scattering theory, ||, |3J], ||, f26|| . The 
condition of quasi-exact solvability imposes an additional constraint on the Lie algebra 
and hence on the type of operators which are allowed. A Lie algebra of first order differ- 
ential operators g will be called quasi- exactly solvable if it possesses a finite-dimensional 
representation space (or module) M C C°° consisting of smooth functions; this means 
that if ip G J\f and J a G g, then J a (ip) G M . A differential operator T is called quasi- 
exactly solvable if it lies in the universal enveloping algebra of a quasi-exactly solvable 
Lie algebra of differential operators. Clearly, the module M is an invariant space for T, 
i.e., T{J\f) C A/", and hence T restricts to a linear matrix operator on N '. We will call the 
eigenvalues and corresponding eigenfunctions for the restriction T\Af algebraic, since they 
can be computed by algebraic methods for matrix eigenvalue problems. (This does not 
mean that these functions are necessarily algebraic in the traditional pure mathematical 
sense.) Note that the number of such algebraic eigenvalues and eigenfunctions equals 
the dimension of M . So far we have not imposed any normalizability conditions on the 
algebraic eigenfunctions, but, if they are normalizable, then the corresponding algebraic 
eigenvalues give part of the point spectrum of the differential operator. 

It is of great interest to know when a given differential operator is in Lie algebraic 
or quasi-exactly solvable form. There is not, as far as we know, any direct test on the 
operator in question that will answer this in general. Consequently, the best approach to 
this problem is to effect a complete classification of such operators under an appropriate 
notion of equivalence. In order to classify Lie algebras of differential operators, and 
hence Lie algebraic and quasi-exactly solvable Schrodinger operators, we need to precisely 
specify the allowable equivalence transformations. 

Definition 1 Two differential operators T(x) and T{x) are equivalent if they can be 
mapped into each other by a combination of change of independent variable, 

x = tp(x), (4) 

and gauge transformation 

T(x) = e a{x) ■T{x)-e- a{x \ (5) 

The transformation (§)-(0) have two key properties. First, they respect the commutator 
between differential operators, and therefore preserve their Lie algebra structure. In par- 
ticular, if g is a finite-dimensional Lie algebra of first-order differential operators then the 
transformed algebra g = {T | T G g} is a Lie algebra isomorphic to g. Moreover, if M is a 
finite-dimensional g-module then A" = {e a ^ f(x)\ x=(p -i(x) \ f G Af} is a finite-dimensional 
g-module. In other words, if g is quasi-exactly solvable so is g. It immediately follows 
that the transformation (§)-(§) preserves the class of quasi-exactly solvable operators: 
in other words, if T(x) is quasi-exactly solvable with respect to g then the transformed 
operator T(x) will be quasi-exactly solvable with respect to the Lie algebra g ~ g. 

Second, preserve the spectral problem (|l[) associated to the differential oper- 

ator T, so that if ip(x) is an eigenfunction of T with eigenvalue A, then the transformed 
(or "gauged") function 

ip(x) = e a ^ip{x), where x = <p(x), (6) 
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is the corresponding eigenfunction of T having the same eigenvalue. Therefore, this notion 
of equivalence is completely adapted to the problem of classifying quasi-exactly solvable 
Schrodinger operators. The gauge factor = e a ^ in (|5|) is not necessarily unimodular, 
i.e., a{x) is not restricted to be purely imaginary, and hence does not necessarily preserve 
the normalizability properties of the associated eigenfunctions. Therefore, the problem 
of normalizability of the resulting algebraic wave functions must be addressed. 

Let us summarize the basic steps that are required in order to obtain a complete 
classification of quasi-exactly solvable operators and their algebraic eigenfunctions: 

1. Classify finite-dimensional Lie algebras of differential operators. 

2. Determine which Lie algebras are quasi-exactly solvable. 

3. Solve the equivalence problem for differential operators. 

4. Determine normalizability conditions. 

5. Solve the associated matrix eigenvalue problem. 

In this survey we will concentrate on the solution of the first three problems, referring 



the reader to [[IT]] for the solution of the normalizability problem in the one- dimensional 

case. 



3 Equivalence of Differential Operators 

Consider a second-order linear differential operator 

defined on an open subset M C M d . We are interested in studying the problem of when 
two such operators are equivalent under the combination of change of variables and gauge 
transformation (f|), ([5]). Of particular importance is the question of when T is equivalent 
to a Schrodinger operator, which we take to mean an operator Ti, = —A + V(x), where A 
denotes either the flat space Laplacian or, more generally, the Laplace-Beltrami operator 
over a curved manifold. This definition of Schrodinger operator excludes the introduction 
of a magnetic field, which, however, can also be handled by these methods. As we shall see, 
there is an essential difference between one-dimensional and higher dimensional spaces in 
the solution to the equivalence problem for second-order differential operators, because in 
higher space dimensions not every second-order differential operator is locally equivalent 
to a Schrodinger operator —A + V(x), where A is the flat space Laplacian. 

Explicit equivalence conditions were first found by E. Cotton, M, in 1900. Since the 
symbol of a linear differential operator is invariant under coordinate transformations, 
we begin by assuming that the operator is elliptic, meaning that the symmetric matrix 
G(x) = (g l: >(x)) determined by the leading coefficients of — T is positive-definite. Owing 
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to the induced transformation rules under the change of variables (|7|), we can interpret 
the inverse matrix G(x) = = (gij(x)) as defining a Riemannian metric 

d 

ds 2 = g i j{x)dx l dx\ (8) 
»J=i 

on the subset M C M. d . We will follow the usual tensor convention of raising and lowering 
indices with respect to the Riemannian metric (|8|). We rewrite the differential operator 
(0) in a more natural coordinate-independent form as 

n 

T=-J2 ^'( V * - A i)Wj - Aj) + V, (9) 

where Vj denotes covariant differentiation using the Levi-Civita connection associated to 
the metric ds 2 . The vector A(x) = (A 1 (x), . . . ,A d (x)) can be thought of as a (general- 
ized) magnetic vector potential; in view of its transformation properties, we define the 
associated magnetic one-form 

d 

A = Y i A i {x)dx i . (10) 

i=i 

(Actually, to qualify as a physical vector potential, A must be purely imaginary and sat- 
isfy the stationary Maxwell equations, but we need not impose this additional physical 
constraint in our definition of the mathematical magnetic one-form (|10"D -) The explicit 
formulas relating the covariant form ([]) to the standard form (I7p of the differential opera- 



tor can be found in [jig] . Each second-order elliptic operator then is uniquely specified by 
a metric, a magnetic one-form, and a potential function V(x). In particular, if the mag- 
netic form vanishes, so A = 0, then T has the form of a Schrodinger operator T = — A+V, 
where A is the Laplace-Beltrami operator associated with the metric (§). 

The application of a gauge transformation (|5]) does not affect the metric or the po- 
tential; however, the magnetic one-form is modified by an exact one-form: A i— > A + da. 
Consequently, the magnetic two-form Q = dA, whose coefficients represent the associated 
magnetic field, is unaffected by gauge transformations. 

Theorem 2 Two elliptic second-order differential operators T and T are (locally) equiv- 
alent under a change of variables x = tp(x) and gauge transformation if and only if 
their metrics, their magnetic two-forms, and their potentials are mapped to each other 

v *(ds 2 ) = ds 2 , (p*(U) = Q, (p*(V) = V. (11) 

(Here if* denotes the standard pull-back action of if on differential forms; in particular, 

<p*(y) = v°<p.) 

In particular, an elliptic second-order differential operator is equivalent to a Schro- 
dinger operator —A + V if and only if its magnetic one-form is closed: dA = Q = 0. 
Moreover, since the curvature tensor associated with the metric is invariant, T will be 
equivalent to a "flat" Schrodinger operator if and only if the metric ds 2 is flat, i.e., has 
vanishing Riemannian curvature tensor, and the magnetic one-form is exact. 
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In the one-dimensional case, every metric is automatically flat and all 1-forms are exact. 
Hence the previous theorem implies that every elliptic second-order differential operator 
is equivalent to a flat Schrodinger operator. See H] for explicit formulas for the change 
of variables, the gauge factor and the potential. 



4 Lie Algebras of Differential Operators 

In this section, we summarize what is known about the classification problem for Lie 
algebras of first order differential operators. Any finite-dimensional Lie algebra g of first 
order differential operators in WL d (with K = E or K = C) has a basis of the form 

J 1 = v 1 + rpix), . . . , J r = v r + r] r (x), 

r +1 = f\x),...,J r+s = f s {x), (12) 

cf. (U). Here v 1 , . . . , v r are linearly independent vector fields spanning an r-dimensional 
Lie algebra f). The functions / 1 (x), . . . , f s (x) define multiplication operators, and span 
an abelian subalgebra JF of the full Lie algebra g. Since the commutator [v l , /•?'] = v* (/■?') 
is a multiplication operator, which must belong to g, we conclude that 1) acts on J 7 , 
which is a finite-dimensional Ij-module (representation space) of smooth functions. The 
functions i] a (x) must satisfy additional constraints in order that the operators (0) span 
a Lie algebra; we find 

[ v * + v \ v J + rf] = [v\ v J ] + vV) - v J '(r/). (13) 

Now, since f) is a Lie algebra, [v l , v J ] = J2k c k vfc > where are the structure constants of 
f). Thus the above commutator will belong to g if and only if 

v'^-v^M-^c^^. (14) 

k 

These conditions can be conveniently re-expressed using the basic theory of Lie algebra 
cohomology, P^| . Define the one-cochain F on the Lie algebra of vector fields f) as the 
linear map F: f) — > C°° = C°°(K a! ) which satisfies (F;v a ) = 7] a . Since we can add in 
any constant coefficient linear combination of the f b, s to the 7/ a 's without changing the 
Lie algebra g, we should interpret the 7/ a 's as lying in the quotient space and 
hence regard F as a C°° /J 7 - valued cochain. In view of (|T3"D, the collection of differential 
operators ([12]) spans a Lie algebra if and only if the cochain F satisfies 



v(F;w) - w(F; v) - (F ; [v,w]) G T, Vv.wefj. (15) 
The left hand side of ([T5| ) is just the evaluation (6xF;v, w) of the coboundary of the 



1-cochain F, hence (15) expresses the fact that the cochain F must be a C°°/jF-valued 
cocycle, i.e., <5i_F = 0. A 1-cocycle is itself a coboundary, written F = S a for some a(x) G 
C 00 , if and only if (F ; v) = v(cr) for all v G f), where v(cr) is considered as an element of 
C°° jT. It can be shown that two cocycles will differ by a coboundary 5q<j if and only if the 
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corresponding Lie algebras are equivalent under the gauge transformation (|5|). Therefore 
two cocycles lying in the same cohomology class in the cohomology space H 1 ^, C°°/JF) = 
Ker 5i/Im 5q will give rise to equivalent Lie algebras of differential operators. In summary, 
then, we have the following fundamental characterization of Lie algebras of first order 
differential operators: 

Theorem 3 There is a one-to-one correspondence between equivalence classes of finite 
dimensional Lie algebras q of first order differential operators on M and equivalence 
classes of triples (I), J 7 , [F]), where 

1. f) is a finite- dimensional Lie algebra of vector fields, 

2. T C C°° is a finite- dimensional \)-module of functions, 

3. [F] is a cohomology class in H l (t), C°° /J 7 ). 

Based on Theorem |3], there are three basic steps required to classify finite dimensional 
Lie algebras of first order differential operators. First, one needs to classify the finite 
dimensional Lie algebras of vector fields f) up to changes of variables; this was done 
by Lie in one and two complex dimensions, [071, and by the authors in two real 



dimensions, |T6| , under the regularity assumption that the Lie algebra has no singularities 
- not all vector fields in the Lie algebra vanish at a common point . Secondly, for each of 
these Lie algebras, one needs to classify all possible finite-dimensional f)-modules T of C°° 
functions. Finally, for each of the modules J 7 , one needs to determine the first cohomology 
space i? 1 (f), C°° /J 7 ). As the tables at the end of |15[ indicate, the cohomology classes are 



parametrized by one or more continuous parameters or, in a few cases, smooth functions. 

It is then a fairly straightforward matter to determine when a given Lie algebra satisfies 
the quasi-exact solvability condition that it admit a non-zero finite-dimensional module 
M C C°°. Remarkably, in all known cases, the cohomology parameters are "quantized", 
the quasi-exact solvability requirement forcing them to assume at most a discrete set of 
distinct values. This intriguing phenomenon of "quantization of cohomology" has been 
geometrically explained in the maximal cases in terms of line bundles on complex surfaces 
in P. 

In one dimension, there is essentially only one quasi-exactly solvable Lie algebra of 
first order differential operators: 

Theorem 4 Every (non-singular) finite- dimensional quasi-exactly solvable Lie algebra of 
first order differential operators in one (real or complex) variable is locally equivalent to 
a subalgebra of one of the Lie algebras 

Q n = Span{ d x , xd x , x 2 d x -nx, 1}~ gl(2), (16) 

where n G N. For g n , the associated module M = V n consists of the polynomials of degree 
at most n. 



Turning to the two-dimensional classification, a number of additional complications present 
themselves. First, as originally shown by Lie, there are many more equivalence classes of 



LIE ALGEBRAS AND PARTIAL INTEGRABILITY 



9 



finite-dimensional Lie algebras of vector fields, of arbitrarily large dimension and some- 
times depending on arbitrary functions as well. Moreover, the classification results in 
M. 2 and C 2 are no longer the same; in fact, the real classification has been completed 



only very recently, |20| . Another complication is that the modules T for the vector field 
Lie algebras are no longer necessarily spanned by monomials, a fact that makes the de- 
termination of the cohomology considerably more difficult. Our classification results for 
finite-dimensional Lie algebras of differential operators in two real or complex variables, 
13fl , fll5| , p0| , are summarized in the Tables appearing, respectively, in p0| and p~4 ] . 



The classification naturally distinguishes between the imprimitive Lie algebras, for which 
there exists an invariant foliation of the manifold, and the primitive algebras, having no 
such foliation. Each of the Lie algebras appearing in the complex classification has an 
obvious real counterpart, obtained by restricting the coordinates to be real. Moreover, 
as explained in ||20[| , in such cases the associated real Lie algebras of differential opera- 



tors and finite-dimensional modules are readily obtained by restriction. As a matter of 
fact, every imprimitive real Lie algebra of vector fields in M? is obtained by this simple 
procedure. In addition, there are precisely five primitive real Lie algebras of vector fields 
in two dimensions that are are not equivalent under a real change of coordinates to any 
of the Lie algebras obtained by straightforward restriction of the complex normal forms. 
The complete list of these additional real forms along with their canonical complexifica- 
tions appear in Table 4 of f20|| . The maximal (complex) algebras, namely Case 1.4 of 
|2(J (sl(2) 0sl(2)), Case 2.3 (*I(3)), and Case 1.11 (g[(2) X R r ) play an important role in 
Turbiner's theory of differential equations in two dimensions with orthogonal polynomial 
solutions, p9[ . Of the additional real normal forms, Case 4.2 (so (3)) was extensively used 



by Shifman and Turbiner, [33], while Case 4.5 (so(3, 1)) is related to recent examples of 
quasi-exactly solvable Hamiltonians in two dimensions constructed by Zaslavskii, (|2 
[53. 



5 Quasi-Exactly Solvable Operators on the Line 

Let us now specialize to problems in one dimension. In view of Theorem |], we let n G N 
be a nonnegative integer, and consider the Lie algebra Q n spanned by the differential 
operators 

which satisfy the standard sl(2) commutation relations. In this section, we shall use z 
instead of x for the "canonical coordinate", retaining x for the physical coordinate in 
which the operator takes Schrodinger form. Since Q n differs from the Lie algebra g n in 
Theorem [| only by the inclusion of constant functions, any Lie-algebraic operator @ 
for the full algebra Q n is automatically a Lie- algebraic operator for the subalgebra g n . 
Therefore, in our analysis of Lie-algebraic differential operators we can, without loss of 
generality, concentrate on the Lie algebra Q n . 
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Using fllTD , it is readily seen that the most general second-order quasi-exactly solvable 
operator in one space dimension can be written in the canonical form 

m d 2 ( n — 1 ,1 d [ n n{n — 1) „] . , 



tfc 2 P 2 ) dz \ 2^ 12 
where -P(^), Q(z), and -R(z) are (general) polynomials of respective degrees 4, 2, 0, whose 



explicit expression in terms of the constants c a b and c a can be found in [14]. Since the 
module Af is the space V n of polynomials of degree at most n, the algebraic eigenfunc- 
tions of flI8|) will, in the z-coordinates, just be polynomials Xk{ z ) £ Pn- In terms of the 
standard basis v^z) = z k , k — 0, . . . , n, the restriction T | V n takes the form of a penta- 
diagonal matrix. In summary, for a normalizable one-dimensional quasi-exactly solvable 
operator there are n + 1 algebraic eigenf unctions which, in the canonical z coordinates, 
are polynomials of degree at most n. 

Specializing the solution to the equivalence problem given by Theorem |^ to the oper- 
ator (|18D , we find (cf. |i~4][ ) that if P(z) > the change of variables required to place the 



operator into physical (Schrodinger) form will, in general, be given by an elliptic integral 

x = ^z)=[ Z -pL= (19) 

the corresponding gauge factor being 

^z) = P(zy^ex P y Z -^d y y (20) 

The potential is given by 

n{n + 2) (PP" - \P' 2 ) + 3(n + 1)(QP' - 2PQ') - 3Q 2 
V(x) = ^ 4 J up R, (21) 

where the right hand side is evaluated at z — if' 1 (x). In the physical coordinate, the 
associated algebraic wave functions will then take the form 

ij(x)= f i(^ 1 (x))-x(v~\x)), (22) 

where x{ z ) is a polynomial of degree at most n. 

The canonical form fll8|) of a quasi-exactly solvable differential operator is not unique, 
since there is a "residual" symmetry group which preserves the Lie algebra Q n . Not 
surprisingly, this group is GL(2, R), which acts on the (projective) line by linear fractional 
(Mobius) transformations 

z^w= az + f3 A=( a detA = A = aS-p-f^0. (23) 



jz + 5 ' V 7 8 

To describe the induced action of the transformations (|23|) on the quasi-exactly solvable 
operators QT8|), we first recall the basic construction of the finite-dimensional irreducible 
rational representations of the general linear group GL(2,R). 
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Definition 5 Let n > 0, i be integers. The irreducible multiplier representation p ni of 
GL(2, M) is defined on the space V n of polynomials of degree at most n by the transfor- 
mation rule 

P(w) i ► = (a5 - /3 7 ) 1 (7* + 5)" P f , P e V n . (24) 

The multiplier representation p n .j has infinitesimal generators given by the differential 
operators (|i~7|) combined with the operator of multiplication by n + 2i representing the 
diagonal subalgebra (center) of g[(2,IR). The action (|23|) induces an automorphism of 
the Lie algebra g n , which is isomorphic to the representation P2-1, and, consequently, 
preserves the class of quasi-exactly solvable operators associated with the algebra jj n . 
Moreover, the corresponding gauge action 

f(z) = ( 7 z + 5) n ■ T{w) ■ ( 7 * + 5)~ n , w = ^4 ( 25 ) 

72 + 

will preserve the space of quasi-exactly solvable operators (|18|). Identifying the operator T 
with the corresponding quartic, quadratic and constant polynomials P, Q, R, we find that 
the action (|25| ) of GL(2,R) on the space of quasi-exactly solvable second-order operators 
is isomorphic to the sum of three irreducible representations, p 4 _ 2 © P2-1 © Po,o! the 
quartic P(z) transforms according to P4-2, the quadratic Q(z) according to P2-1, while 
R is constant. Finally, the associated module, which is just the space of polynomials V n , 
transforms according to the representation p n - 

Using the action of GL(2,R), we can place the gauged operator (0) into a simpler 
canonical form, based on the invariant-theoretic classification of canonical forms for real 



quartic polynomials, [Q, [fP7 | . 



Theorem 6 Under the representation P4-2 o/GL(2, M), every nonzero real quartic poly- 
nomial P(z) is equivalent to one of the following canonical f orms: 

~z 2 )(l -k 2 z 2 ), v{l- z 2 ){l-k 2 + h 2 z 2 ), u{l + z 2 )(l + h 2 z 2 ) } 

v(z 2 -l), v(z 2 + l), uz 2 , u(z 2 + l) 2 , z, 1, (26) 

where v and < k < 1 are real constants. 



The solution to the normalizability problem, which the interested reader can find in fL7 
begins with a detailed analysis of the elliptic integral (fL9f). It is found, first of all, that 
the class of quasi-exactly solvable potentials naturally splits into two subclasses — the 
periodic potentials, which are never normalizable, and the non-periodic potentials, which 
are sometimes normalizable. Tedious but direct calculations based on (|i~9D, (pop, (pT]), 
and (p2|), produce then the explicit change of variables, the potential, and the form of the 
algebraic eigenfunctions for the above normal forms in physical coordinates. Each of the 
classes of potentials is a linear combination of four elementary and/or elliptic functions, 
plus a constant which we absorb into the eigenvalue. The potentials naturally divide into 
two classes, which are listed in the following two Tables. In each case, the four coefficients 
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are not arbitrary, but satisfy a single complicated algebraic equation and one or more 
inequalities; see |17[ for the details. 

First, the periodic quasi-exactly solvable potentials correspond to the cases when the 
real roots of P (if any) are simple, of which there are five cases in (26). The explicit 
formulas for the corresponding potentials follow; they depend on two real parameters 
a > and k e (0, 1). In the first three cases, the corresponding potentials are written in 
terms of the standard Jacobi elliptic functions of modulus k, |S], vol. 2]. Also, as remarked 
above, the coefficients A, B, C, D are not arbitrary, although the explicit constraints are 
too complicated to write here. 



Periodic Quasi-Exactly Solvable Potentials 

1. dn -2 a x(A snax + B) + cn~ 2 a x(C snax + D) 

2. dn -2 a x(A cnai + B) + sn~ 2 a x(C cn a x + D) 

3. Acnaxsnax + B cn 2 a x + C duT 2 a x(cn axsnax + D cn 2 a x) 

4. A sin 2 ax + B sin ax + C tan a x sec ax + D sec 2 a x 

5. A cos 4a x + B cos la x + C sin 2a x + D sin 4a x 

Note that the potentials in cases 1, 2 and 4 have singularities unless C = D = 0. In 
Cases 3 and 5, the potential has no singularities, reflecting the fact that in these cases 
P(z) has no real roots. Case 3 includes the Lame equation, ||, vol. 3], and Case 4 with 
A = B = is the trigonometric Scarf potential. 

The non-periodic potentials correspond to the cases with one or two multiple roots. 
The explicit formulas for the corresponding potentials follow (again, a > is a real 
constant). 



Non-periodic Quasi-Exactly Solvable Potentials 

1 . A sinh 2 ax + B sinh a x + C tanh a x sech ax + D sech 2 a x 

2. A cosh 2 ax + B cosh ax + C coth a x csch ax + D csch 2 a x 

3. Ae 2ax + Be ax + Ce~ ax + De~ 2ax 



4. Ax b + Bx* + Cx 2 + 



D 



5. Ax 4 + Bx 3 + Cx 2 + Dx 



In cases 2 and 4, the potential has a singularity at x = unless C + D = (Case 2) 
or D = (Case 4). The nonsingular potentials in Case 4 are the anharmonic oscillator 
potentials discussed in detail in |34j], fl32"| . The algebraic constraints satisfied by the 
coefficients are given in [T7 . 



According to Turbiner, 0], a potential is exactly solvable if it does not explicitly 
depend on the discrete "spin" parameter n, since, in this case, one can find representation 
spaces of arbitrarily large dimension and thereby (if the algebraic eigenfunctions are 
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normalizable) produce infinitely many eigenvalues by algebraic methods. Note that since 
the gauge transformation (|19D , (p0|) , can explicitly depend on n, exact solvability cannot 
be detected in the canonical coordinates, but depends on the final physical form of the 
operator. The exactly solvable nonperiodic potentials are characterized by the condition 
A = B = 0, and, in case 3, C = D = 0. In Case 2, there is an additional inequality, 
\C\ < D + \oi 2 1 to be satisfied. The exactly solvable potentials include the (restricted) 
Poschl-Teller and Scarf potentials (Case 1), the Rosen-Morse II potential (Case 2), the 
Morse potential (Case 3), the radial harmonic oscillator (provided D — 1(1 + 1), I e N) 
(Case 4), and the harmonic oscillator (Case 5). 

Analysis of the explicit formulas for the eigenfunctions based on ( p2]) yields a complete 
set of conditions for the normalizability of the non-periodic potentials, which can be 
written explicitly in terms of the coefficients of the quadratic polynomial Q(z). It is 
then possible to deduce explicit, general normalizability conditions on the Lie-algebraic 
coefficients c a b and c a by using the fact that such conditions must be invariant under the 
action of the group GL(2, R). Therefore, normalizability conditions can be written in 
terms of the classical joint invariants and covariants of the pair of polynomials P, Q. See 
|T7| for the details and a complete list of invariant normalizability conditions. 



6 Two- Dimensional Problems 

There are a number of additional difficulties in the two-dimensional problem which do not 
appear in the scalar case. First, there are several different classes of quasi-exactly solvable 
Lie algebras available. Even more important is the fact that, according to Theorem ||], 
there are nontrivial closure conditions which must be satisfied in order that the magnetic 
one-form associated with a given second-order differential operator be closed, and hence 
the operator be equivalent under a gauge transformation (|5|) to a Schrodinger operator. 
Unfortunately, in all but trivial cases, the closure conditions associated with a quasi- 
exactly solvable operator ([|) corresponding to the generators of one of the quasi-exactly 
solvable Lie algebras on our list are nonlinear algebraic equations in the coefficients c a b, c a , 
Co, and it appears to be impossible to determine their general solution. Nevertheless, there 
are useful simplifications of the general closure conditions which can be effectively used 
to generate large classes of planar quasi-exactly solvable and exactly solvable Schrodinger 
operators, both for flat space as well as curved metrics. 

Suppose that the Lie algebra g is spanned by linearly independent first-order differ- 
ential operators as in (0). The closure conditions dA = for a Lie-algebraic second-order 
differential operator (Q) are then equivalent to the solvability of the system of partial 
differential equations 

r 

2^2c ab r] b + c a , 

. 6=1 

i = l,...,d, (27) 

for a scalar function r(x), given by r = 2a + \ logdet^-) in terms of the gauge factor 
e u required to place the operator in Schrodinger form. The closure conditions ( p7[) are 
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extremely complicated to solve in full generality, but a useful subclass of solutions can 
be obtained from the simplified closure conditions 



2r] a = k a , 



r. 



(28) 



where k l , . . . , k r are constants. Any solution t(x) of equations ( p8|) will generate an infin- 
ity of solutions to the full closure conditions (0), with c a b arbitrary, and c a = J2b c abk b - 
The case k a = and g semi-simple was investigated in f29fl . Although the simplified 
closure conditions can be explicitly solved for such Lie algebras, with the exception of 
so (3), their solutions are found to generate quasi-exactly solvable Schrodinger operators 
that are not normalizable, and hence of limited use. Note that even when the simplified 
closure conditions do not have any acceptable solutions, the full closure conditions 
may be compatible and may give rise to normalizable operators, as shown in [|18[ and 
2"0fl . Here we shall limit ourselves to a few examples. 

Consider, in the first place, the Lie algebra g ~ sl(2)©5[(2) spanned by the first-order 
differential operators 



J 1 



J 1 = d x , J 2 = 0, 
yd y , J 5 = x 2 d x 



yi 

nx, 



J 3 - 
J 6 



xO x , 

= y 2 9 y 



my, 



(29) 



where n, m G N. The particular choice 



(Cab 



( 2 


i 











1\ 


1 


2 








1 











3 




















3 











1 








1 


1 


Vi 











1 


1/ 



0,0,-(l + 4n),-(l + 4m),0,0 



c 



+ m + n 



(30) 



(31) 



of Lie-algebraic coefficients lead to a quasi-exactly solvable Hamiltonian with Riemannian 
metric 



.ii 



9 



22 



l+y 2 )(2 + y 2 



g - = (1 + x *)(2 + x 2 ), g 12 = (1 + x 2 )(l + y 2 ) 
which has complicated curvature, and potential 

(l + 2n)(3 + 2n) (l + 2m)(3 + 2m) 



(32) 



AV 



-y 



1 + x 2 1 + y 2 

17 + 12y 2 -y 4 + 2xy(6 + 5y 2 ) 



3 + x 2 



y 



+ 



5(3 + 2xy)(l + y 2 )(2 + y 2 
(3 + x 2 + y 2 ) 2 



(33) 
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Consider next the so (3, 1) Lie algebra generated by the operators 



J 1 = d x , J 2 = d y , J 3 = xd x + yd y , J 4 = yd x - xd, 



J 5 = (x 2 — y 2 )d x + 2xyd y — 2nx, J 6 = 2xyd x + (y 2 — x 2 )d y — 2ny, 



(34) 



with n G N a positive integer. Notice that this algebra, although of course inequivalent 
to any of the real normal forms in the classification of p(J , is comp/ex-equivalent to the 
sl(2) © sl(2) Lie algebra in the previous example under the complex change of variables 
z = x + iy, w = x — iy. 

The second-order differential operator 



— T 



a 



t4\2 



+ A 



(ry+ (j b y +Ap(i + 2 n y -4 7 



(35) 



satisfies the closure conditions, and is therefore equivalent to a Schrodinger operator on 
the manifold with contravariant metric tensor {g 1 ^) given by 



9 
9 
9" 



ii 



12 



. 22 



a + (3 x 2 + 7£/ 2 + A(x 2 + y 2 ) 2 , 
a + -fx 2 + py 2 + \(x 2 + y 2 ) 2 . 



The Gaussian curvature is 

(-P + 3 7 )(a 2 + A 2 r 8 ) + 2(/3 7 + 4aA)r 2 (a + Ar 4 ) + 2a\(5(3 + 7 )r 4 



(36) 



(37) 



(a + 7 r 2 + Ar 4 ) 2 

with r 2 = x 2 + y 2 . If the parameters a, (3, 7 and A are positive, then the metric is 
non-singular and positive definite for (x,y) ranging over all of M. 2 . The fact that the 
closure conditions are satisfied guarantees the existence of a gauge factor fi such that 
H = \iT fi^ 1 is a Schrodinger operator; see |20[ for an explicit formula for the gauge 
factor, which is too complicated to be given here. The expression for the potential V is 



AV 



16 a 13 n (1 + n) + r 2 [f3 2 (3 + 16 n + 16 n 2 ) - 4 a A (3 + 8 n + 4 n 2 )] 



+ 



a + (3r 2 + Ar 4 
5Q3-7) (4 a A - 7 2 ) + 3A(2/3 7 - 3 7 2 + 4aA)r 2 

A(a + 7 r 2 + Ar 4 ) 
5Q3-7) (4a A - 7 2 ) (a + 7 r 2 ) 
A(a + 7 r 2 + Ar 4 ) 2 ' 



(3? 



Since the potential is a function of r only, it is natural to look for eigenf unctions of H 
which depend on r only. When this is done, it can be shown that one ends up with an 
effective Hamiltonian on the line which is an element of the enveloping algebra of the 
standard realization of sl(2, M) in one dimension. Thus, no new quasi-exactly solvable 
one-dimensional potentials are obtained by reduction of the above quasi-exactly solvable 
so(3, 1) potential. This lends additional support to the observation that reduction of 
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two-dimensional quasi-exactly solvable Schrodinger operators does not lead to any new 
one-dimensional quasi-exactly solvable operators. 

Non semi-simple Lie algebras can also yield interesting examples of quasi-exactly 
solvable potentials. Indeed, let g be the Lie algebra spanned by the first-order differential 
operators 

J 1 = d x , J 2 = d y , J 3 = xd x , J 4 = xd y , J 5 = yd y , (39) 

and 

J 6 = x 2 d x + (r - l)xyd y - nx, J 6+i = x i+1 d y , i = 1, . . . , r - 2. (40) 

The module TV is now spanned by the monomials x l yi with i + (r — l)j < n. For m G N, 
a, j3 > 0, the Schrodinger operator with metric 

g n = ax 2 + (3, g 12 = (1 + m)axy, g 22 = (ax 2 + (3) m + a(l + m) V, (41) 

and potential 

Xa(3(l + m) 2 (ax 2 + (3) m ^ 1 . n/ lN //10 . 

(ax 2 + p) i+m + ap(l + myy 1 

is normalizable and quasi-exactly solvable with respect to q, provided that the parameter 
A is large enough. The metric in this case has constant negative Gaussian curvature k = 
—a. Furthermore, since the potential V does not depend on the cohomology parameter 
n, the above Hamiltonian is exactly solvable. Moreover, the potential is also independent 
of r; hence we have constructed a single exactly solvable Hamiltonian which is associated 
to an infinite number of inequivalent Lie algebras of arbitrarily large dimension. 



7 Matrix Schrodinger operators 

We shall outline in this Section how to extend the notion of quasi-exact solvability to 
matrix Hamiltonians in one dimension. By definition, a matrix Schrodinger operator (or 
matrix Hamiltonian) is a iV x N matrix of second-order differential operators of the form 
7i = —d 2 + V(x), where V(x) is a. N x N Hermitian matrix of functions. If we drop the 
restriction that V be Hermitian, then 7i will be called a Schrddinger-like operator. Here 
we shall be mainly interested in spin 1/2 matrix Hamiltonians, that is in the case N = 2; 
see JTDJ for the case of arbitrary N. 

The notion of equivalence we shall use for matrix differential operators is the same as 
in the scalar case (cf. Definition [I]), where the gauge factor e "^ is replaced now by an 
arbitrary invertible matrix U(x). Suppose, as in the previous sections, that a Schrodinger 
operator 

H(x) = U(z)-T(z)-U(z)-\ x = <p{z), (43) 

is equivalent to a matrix differential operator T(z) preserving a finite-dimensional sub- 
space M C C°°(R) © C°°(M) of smooth two-component wave functions, i.e., T(M) C 
M . It follows that TC will restrict to a linear operator in the finite-dimensional space 
U ■ Af\ z=v -i( x ), and therefore dimA/" eigenfunctions of Tt will be computable by purely 
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algebraic methods. We shall say that Ti is quasi- exactly solvable; as in the scalar case, 
the algebraic eigenfunctions of 7i should obey appropriate boundary conditions, like for 
instance square integrability. 

Thus, to find examples of quasi-exactly solvable matrix Hamiltonians we have to solve 
the following two problems: 

1. Characterize all second-order matrix differential operators leaving invariant a finite- 
dimensional space of smooth functions 

2. Find when a second-order matrix differential operator satisfying the previous con- 
dition is equivalent to a matrix Schrodinger operator 

The first of these problems is clearly too general. In the scalar case, it is replaced by 
the simpler (but still highly nontrivial) problem of constructing second-order differential 
operators belonging to the enveloping algebra of a Lie algebra of quasi-exactly solvable 
first order differential operators, which in turn leads to the classification of such Lie 
algebras. In the matrix case this is still too general, since Lie superalgebras and even 
more general algebraic structures, whose complete classification is not available yet, come 
into play, 0], §, [0. 



To simplify matters, the idea is to fix a concrete (but still general enough) finite- 
dimensional subspace of smooth wave functions and classify all second-order matrix dif- 
ferential operators leaving this subspace invariant. Following Turbiner, [38fl, and Brihaye 
et al, [|1|, 0, we choose as our invariant subspace the space of all wave functions with 
polynomial components. More precisely, let 

be the vector space of wave functions ^(z) = (ipi(z), ^(z)) , with ipi(z) G V m and 
i>2(z) G V n polynomials of degrees less than or equal to m and n, respectively. We shall 
say that a matrix differential operator T is a polynomial vector space preserving (PVSP) 
operator if it preserves the vector space V m ,n for some non-negative integers m and n. 
We shall also denote by V$ n the vector space of all PVSP matrix differential operators of 
differential order less than or equal to k. The problem is then to characterize all second- 
order PVSP operators, or in other words to describe the space V$ n in as concise a way 
as possible. 

To this end, let A = n — m > 0, and consider the 2x2 differential operators 

T+ = diag(J+_ A , J+), T° = diag(J°_ A , J° n ), 
T~ = diag(J~, J"), J = \ diag(n + A, n); 
Q a = z a a-, Q a = q a {n, A) a+ , « = 0,...,A, (45) 



withQ 



q a (n,A)=l[(zd z -n + A-k)d: and a + = (a'f = ( ^ ^ ] . ( 10) 

k=l 




^Here, and in what follows, we have adopted the convention that a product with its lower limit greater 
than the upper one is automatically 1. 
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Note that, in contrast with the scalar case, for A > 1 the operators Q a have differential 
order grater than one. 

It can be easily checked that the 6 + 2A operators in ( [45]) preserve V n -A,n, as does 
any any polynomial in the above operators. We now introduce a ^-grading in the set 
T>2x2 of 2 x 2 matrix differential operators as follows: an operator 



T= " " , (47) 




where a, b, c and d are differential operators, is even if b = c = 0, and odd if a = d = 0. 
Therefore, the T's and J are even and the Q's and Q's odd. This grading, combined 
with the usual product (composition) of operators, endows V 2x2 with an associative 
superalgebra structure. We can also construct a Lie superalgebra in V 2x2 by defining a 
generalized Lie product in the usual way: 

[A, B] s = AB — (-l) degAdegB 5A (48) 

However, this product does not close within the vector space spanned by the operators in 



(fig), except for A = 0, 1 (see [jnj for the explicit commutation relations). More precisely, 
for A = 1 the underlying algebraic structure is the classical simple Lie superalgebra 
osp(2,2), [|31|, [j38] , whereas for A = it is P)i © s[(2), where f)i is the 3-dimensional 



Heisenberg superalgebra. As remarked in M, in the latter case we can alternatively 
leave the grading aside and replace J = 1 by J = 03, ending up with the Lie algebra 
sl(2) © sl(2). For all other values of A, the Lie superalgebra of V 2x2 generated by the 
operators fl45|) is infinite-dimensional, |L0] . 



In spite of the fact that in general they don't generate a finite-dimensional superalge- 
bra, the operators (|45"D are the building blocks in the construction of second-order matrix 
differential operators preserving V m ,n] in fact, any such operator must be a polynomial 



in the operators ([45]), a fact first stated by Turbiner, p8| : 



Theorem 7 Let A = n — m > 0, m > 2, and let T be a second-order 2x2 matrix 
differential operator leaving the space V m>n invariant. Then 

T = p 2 (T*) + Jp 2 (T*) + ]T Q a p a 2 _ A (T e ) + E Q a p 2 (T £ ), (49) 

with pk, pk, Pk> an d Pt polynomials of degree less than or equal to k. If A = 0, the 
above formula is still valid with J replaced by J = 0-3, while if A > 2 every p% must be 
identically zero. 



See [ID for a proof, as well as for the generalization of this result to the case of N x iV 
matrix differential operators of arbitrary order. 

Turning now to the second of the two fundamental problems listed at the begininng 



of this section, we have the following fundamental Theorem, [JI0 |: 
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Theorem 8 Let T be a PVSP operator in V$ n , with n > m > 2. Then T is equivalent 
to a Schrddinger-like operator if and only if it is of the form 

- T = Pi (z) d 2 + A^z) d z + A (z) , (50) 

with Pa{z) a polynomial of degree at most 4- The operator T is equivalent to a Schrodinger 
operator —d 2 + V(x) if and only if ([57Jj holds, and in addition we have: 

1. The physical coordinate x = <p(z) is given by 

r z ds , , 

x = I -_; (51) 

/p 4 {s) 



2. There is an invertible matrix U (x) satisfying the differential equation 

U x = UA, (52) 

with A given by^\ A(x) = - ~p' 4 )\ z=v -i {x) ; 

3. The matrix 

V = UWIJ- 1 , (53) 
with W = —Aq oip~ l + A 2 + A x , is Hermitian. 

The eigenfunctions of the Hamiltonian H are of the form \I/(x) = U(x)<& (ip~ l (x)), where 
is an eigenfunction ofT and ^ must satisfy appropriate boundary conditions. 

Using the previous theorem, it can be shown that a matrix Schrodinger operator equiv- 
alent to a PVSP operator T G V^Ia n with A > 1 can be transformed by a constant 
gauge transformation into a Schrodinger operator whose potential matrix is triangular. 
Since triangular potentials correspond to essentially uncoupled systems, we can restrict 
ourselves without loss of generality to the cases A = and A = 1. 

The above theorems tell us, in principle, how to find examples of quasi-exactly solv- 
able matrix Hamiltonians possessing a finite number of algebraically computable wave 
functions. The steps in the construction can be summarized as follows: 

1. Construct the most general PVSP operator of the form (fl9|). This operator will 
depend on a finite number of parameters, namely the coefficients of the arbitrary 



polynomials appearing in ( f49|) 



2. Compute the indefinite integral ( 5"I ) to express the physical coordinate x in terms 
of the "canonical coordinate" z 

3. Solve the matrix linear differential equation ( [52]) to compute the gauge factor U(x) 



^Here and in what follows, a subscripted x denotes derivation with respect to x, while derivatives 
with respect to z will be denoted with a prime '. 
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4. Compute the potential matrix V(x) using equation (|53|) 

5. Choose the values of the free parameters in such a way that the matrix V(x) com- 
puted in the previous step is Hermitian, and that the phyisical wave functions 
satisfy the appropriate boundary conditions 



Of the previous five steps, only the first two present practical problems. Indeed, the first 
step requires the computation and inversion of an elliptic integral, while the second one 
involves the solution of a linear matrix differential equation with variable coefficients. 

The first of the latter two problems can be solved, as in the scalar case, by noting 
that there is still a "residual" GL(2, M) symmetry which can be used to put p^ into one 
of the canonical forms (|2~6"D, for which the elliptic integral ( j51"D reduces to an elementary 
function. We refer the interested reader to |ITJ for the complete story. As to the second 
problem (the solution of the matrix differential equation (0)), the strategy is to impose 
additional conditions on the matrix A that enable the explicit integration of (|52"|), without 
being overly restrictive. 

This can be done as follows. First of all, let us rewrite ([52]) in the z variable as 



where 



U' = uA, (54) 



A = AO0) = 1 ( ^ _ 1^ = . (55) 



a/P4 2p 4 2 

Suppose that A satisfies 



[A(z), / A{s)ds] = 0, (56) 



for some zq G M. If this condition holds, we shall say that we are in the commuting case. 
(It can be shown that ( |56|) is indeed verified by the quasi-exactly solvable Schrodinger 
operator found by Shifman and Turbiner, [53].) In the commuting case, we can readily 



integrate the gauge equation (p4|), obtaining the following general solution: 

U(z) = U exp [ Z A{s)ds, where U E GL(2,C). (57) 

It can be easily shown that the most general 2x2 matrix M(z) satisfying fl56|) must be 
of the form 

M(z) = f(z)M + g(z), (58) 

with Mq a constant matrix and /, g scalar functions. Using this simple observation and 
the explicit formulas for the operators (|45|) for the cases A = 0,1, one readily computes 
the matrix A(z) corresponding to each of the canonical forms of p&{z). It turns out that 
the first three canonical forms (|26| ) of lead to uncoupled Hamiltonians, and therefore 
can be discarded without loss of generality; see JTIJ] for the explicit formulas for the 
remaining canonical forms. 
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In the non-commuting case (that is, when [A, J 2 A] 7^ 0), we may still be able to 
integrate the gauge equation (l54f) explicitly by imposing other constraints on A, as e.g., 
assuming it is uncoupled. Unfortunately, however, no interesting examples of quasi- 
exactly solvable Hamiltonians have been found so far in the non-commuting case. 

Once a solution of the gauge equation flS^f ) has been found, the potential matrix V(x) 
can be immediately computed using fl5"3"|). There are still two conditions that we must 
impose: first, that the potential matrix V(x) in flo^ ) be Hermitian, and secondly, that the 
algebraic wave functions satisfy suitable boundary conditions, which in the examples that 
follow simply reduce to square integrability. These conditions (particularly the Hermitic- 
ity of V) are still too complicated to be solved in full generality, even in the commuting 
case. This is completely analogous to what happened in the two-dimensional scalar case 
(Section |]), where the closure conditions cannot be solved in general, and we are there- 
fore restricted to obtaining particular solutions. In the same spirit, we shall exhibit now 
two examples of quasi-exactly solvable spin 1/2 matrix Hamiltonians, all of which belong 
to the commuting case described above; see |IJ for additional examples and further de- 



tails. Interestingly, both of these examples correspond to the case A = 1; indeed, in the 
commuting case with A = every potential we have obtained turned out to be either 
non-normalizable, singular or diagonalizable by a constant gauge transformation. 
For our first example we consider the PVSP operator 

— T = (T ) 2 + 2a 2 T + + 2(n + 1)T° - 2 JT° + 2&i Q 

- 2&x Q T° - (3n + 1)&! Q - Aa 2 b l Q 1 - (2d + n + ~)J, (59) 

with all the parameters real. In this example = z 2 , so that we are in case 6 of 
Theorem |6|. Solving equation ( 51 ) for z we immediately obtain z = e x . The gauge factor 
reads: 

U(z) = Vz~e a2Z ( C ° SU smu ] : where u = b 1 \ogz. (60) 
v ' V — sin u cos u J 

The potential is given by 

Vj = -(2n + l)a 2 e x + a\ e 2x + (-l) i (a(x) cos(26ix) - (3{x) sin(26!x)), 
v = a(x) sin(26ix) + f3(x) cos(26ix) , (61) 

where j = 1,2, and 

a (x) = -y + a 2 e x , (3{x) = (2n + 1) b x + 2a 2 &i e x . (62) 
It is easily verified that the expected value of the potential is bounded from below, i.e., 
(V,VV) > c\\V\\ 2 , with $6L 2 (M)®L 2 (R), (63) 

for some c£l. (Note, however, that even in this case the amplitude of the oscillations 
of v(x) tends to infinity as x — > +00.) Finally, the condition a 2 < is necessary and 
sufficient to ensure the square integrability of the eigenf unctions \&(x). 
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As our last example, we consider: 

— T = T~T° + 2a x T° + (2a +n--)T~- JT~ + 2b x Q - 2b x Q T° 

- 61 (4a + 3n + 1) g - 4ai&iQi + 2(2a — ai) J, (64) 

where all the coefficients are real, and b\ 7^ 0. Since P4 = z (case 7), we have z = x 2 /4. 
The gauge factor is chosen as follows: 

U(z) = COS \ Z ^ ) . 

V — sin 012: cosoiz / 

The entries of the potential V(x) are given by 



2a (2a -l) 1, 2 , 2 x 2 , / ^,7. , v . M\ 
^2 h 4 ~ & i) x + ( a o cos -g a ( x ) sm -g - ) ' 

b\X 2 b\x 2 
do sin — h at(x) cos — — , (65) 



with j = 1, 2, and a(x) is defined by 

a(x) = — (4ao + 4n + 1 + a±x 2 ) . 

We first note that the potential is singular at the origin unless ao = 0, 1/2. Let us 
introduce the parameter A = 2ao — 1, in terms of which the coefficient of x~ 2 in Vj is 
A(A + 1). If A is a non-negative integer I, we may regard 

(-d 2 x + V(x) -E)$(x) =0, 0<x<oo, (66) 

as the radial equation obtained after separating variables in the three-dimensional Schrod- 
inger equation with a spherically symmetric Hamiltonian given by 

H = -A + U(r) , with U(x) = V(x) - l Jt^R , 

x l 

where A denotes the usual flat Laplacian. Given a non-negative integer I and a spherical 
harmonic Yi m (9,4>), if $ is an eigenfunction for the equation (|66|) satisfying 

lim $(x) =0, (67) 

then 

*(r,M) = — *im(M) 

r 

will be an eigenfunction for H with angular momentum I. If A is not a non-negative 
integer, we shall consider (|66"D as the radial equation for the singular potential U (r) = V(r) 
at zero angular momentum. The potential U(r) is physically meaningful, in the sense 
that the Hamiltonian H admits self-adjoint extensions and its spectrum is bounded from 
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below, whenever A 7^ —1/2, [|TT|, []T7|]. The boundary condition ( pTj) must be satisfied in 
the singular case for all values of A. This boundary condition is verified if and only if 
ao > 0. The expected value of the potential is bounded from below, that is, equation ( |63D 
holds, if and only if 



a 1 

h 



> 1 + V2. 



Finally, the conditions 

a > , ai < , 
guarantee the square integrability of the eigenfunctions 
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